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Traditionally, satellite instruments that measure Earth-reflected solar radiation in the visible and near infrared 
wavelength regions have been calibrated for radiance responsivity in a two-step method. In the first step, the 
relative spectral response (RSR) of the instrument is determined using a nearly monochromatic light source such 
as a lamp-illuminated monochromator. These sources do not typically fill the field-of-view of the instrument nor 
act as calibrated sources of light. Consequently, they only provide a relative (not absolute) spectral response for 
the instrument. In the second step, the instrument views a calibrated source of broadband light, such as a lamp- 
illuminated integrating sphere. The RSR and the sphere absolute spectral radiance are combined to determine the 
absolute spectral radiance responsivity (ASR) of the instrument. More recently, a full-aperture absolute 
calibration approach using widely tunable monochromatic lasers has been developed. Using these sources, the 
ASR of an instrument can be determined in a single step on a wavelength-by-wavelength basis. From these 
monochromatic ASRs, the responses of the instrument bands to broadband radiance sources can be calculated 
directly, eliminating the need for calibrated broadband light sources such as lamp-illuminated integrating spheres. 
In this work, the traditional broadband source-based calibration of the Suomi National Preparatory Project (SNPP) 
Visible Infrared Imaging Radiometer Suite (VIIRS) sensor is compared with the laser-based calibration of the 
sensor. Finally, the impact of the new full-aperture laser-based calibration approach on the on-orbit performance 
of the sensor is considered. 


OCIS codes: (120.0120) Instrumentation, measurement and metrology; (120.028) Remote sensing and sensors; (120.4800) Optical standards and 


testing; (120.5630) Radiometry. 
http://dx.doi.org/10.1364/A0.99.099999 


1. Introduction 


SNPP VIIRS and the successor Joint Polar Satellite System (JPSS) 
VIIRS instruments are designed to provide data continuity with the 
Moderate-resolution Imaging Spectroradiometer (MODIS) 
instruments from the National Aeronautics and Space 
Administration's (NASA’s) Earth Observing System (EOS) [1]. VIIRS 
incorporates a substantial hardware heritage from MODIS, including 
its calibration subsystems [1, 2], which are flight-proven designs. The 


VIIRS sensor is a 22-band instrument with 3 focal planes, one in the 
Visible-Near Infrared (VisNIR), one in the Short-Wave Infrared (SWIR) 
and Mid-Wave Infrared (MWIR), and one in the Long-Wave Infrared 
(LWIR). The wavelengths for the spectral bands extend from 412 nm 
to 12.013 um [3]. In this work, we only consider the VisNIR focal 
plane. There are 9 spectral bands on the VisNIR focal plane, 7 
Moderate spatial resolution "M" bands and 2 higher resolution 
Imaging bands (‘T’), as well as a Day-Night Band (DNB) [3]. There are 
16 filtered Si PIN photodiode detectors in each moderate resolution 


band and 32 Si PIN detectors in each imaging band oriented along the 
track direction. Figure 1 is a schematic diagram (not to scale) of the 
VisNIR focal plane showing the locations of the 7 M-bands, the 2 I- 
bands and the DNB [4]. Nominal center wavelengths and bandpasses 
for the 7 visible (Vis) and near-infrared (NIR) bands discussed in this 
work are given in Table 1. 
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Figure 1. Schematic diagram of the SNPP VIIRS focal plane. 


Table1. The SNPP VIRS Moderate Resolution 
Bands’ nominal center wavelengths and bandwidths 
Center 


Wavelength Bandwidth 





The pre-launch calibration of SNPP VIIRS by the instrument vendor 
closely mimicked the pre-launch calibration of the MODIS 
instruments. A Spectral Measurement Assembly (SpMA), based on a 
lamp-illuminated 0.25 m double grating monochromator, was used to 
measure the VIIRS bands’ RSRs. A lamp-illuminated spherical 
integrating sphere source (SIS), the SIS100 [5], was used to determine 
band gains, converting RSRs into ASRs. The SIS100 radiances, 
combined with the VIIRS sensor bands’ pre-thermal vacuum (pre- 
TVAC) RSRs, provide calibrated band-averaged radiance 
responsivities for the visible and near infrared bands. 

The past decade has seen the development ofa laser-based 
radiometric calibration facility at the National Institute of Standards 
and Technology (NIST), the Facility for Spectral Irradiance and 
Radiance Responsivity Calibrations using Uniform Sources 
(SIRCUS)[6]. SIRCUS incorporates wavelength-tunable lasers coupled 
into sources, typically integrating spheres or collimators, for radiance 
and irradiance responsivity calibrations of instruments. The 
radiometric scale is maintained using calibrated transfer standard 
detectors traceable to primary national radiometric standards at the 
NIST. A travelling version of the NIST SIRCUS facility, T-SIRCUS, was 
used in this study. In this paper, the acronyms SIRCUS and T-SIRCUS 
are used interchangeably in text and figures to represent the travelling 
SIRCUS. 

There are four principal advantages to using the laser-based source 
approach over the historical SoMA/lamp-illuminated integrating 
sphere approach. The first advantage is available flux; the tunable 
laser systems provide tens to hundreds of mW ata particular 
wavelength, three or more orders of magnitude greater flux than 


lamp-monochromator systems. The greater flux enables Out-of-Band 
(OOB) response to be determined over a large dynamic range, 
facilitates scattered light and optical cross-talk characterizations of 
systems, and offers the possibility of solar-level illumination for 
sensors like VIIRS that use the sun as an on-orbit calibration source. 

With the lasers introduced into a SIS, the sensor’s entrance pupil 
can be uniformly illuminated with radiance levels approaching or 
exceeding predicted on-orbit radiances. With the conventional lamp- 
monochromator approach, the radiant flux is too low for full aperture 
illumination. Consequently, conventional sources used to determine 
the sensor’s RSR typically under-fill the instrument’s field-of-view and 
a piece-parts approach must be used to determine the system-level 
RSR. Figure 2 shows the T-SIRCUS and the SpMA illumination 
patterns on the VisNIR focal plane. The SpMA illumination is an image 
of the spectrometer exit slit and fills one detector band ata time on the 
focal plane. The SIRCUS illumination pattern, on the other hand, 
overfills the entire focal plane. Because the output of an integrating 
sphere is imaged at the focal plane, the illumination pattern is very 
uniform, typically to within a couple tenths of a percent over the entire 
illumination pattern. 
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Figure 2. Schematic diagram of the SNPP VIIRS focal plane. 


Second, the output of the T-SIRCUS illuminated SIS viewed by the 
satellite instrument under test is essentially unpolarized compared to 
the highly polarized output from the SpMA. Therefore, light from the 
T-SIRCUS illuminated SIS more closely represents the largely 
unpolarized scattering case experienced by on-orbit remote sensing 
instruments. Third, the radiance of the T-SIRCUS SIS at each 
wavelength is known with very low uncertainties, a factor of 5 or 
more less than the uncertainty in the broad-band lamp-illuminated 
SIS radiance. Finally, the wavelength of the laser sources is known to 
better than 0.01 nm, an order of magnitude better than the 
wavelength uncertainty typical in monochromator-based systems. 
Together, these advantages offer the possibility of calibrating satellite 
sensors with lower uncertainty than is possible using historical 
calibration approaches. 

To evaluate the potential of the new calibration approach for 
satellite sensors, the SNPP Project Management Office at NASA’s 
Goddard Space Flight Center (GSFC) funded a Special Test of VIIRS 
using T-SIRCUS at the SNPP spacecraft manufacturer's facility (Ball 
Aerospace and Technologies Corporation, BATC) in early 2010. This 
test was performed immediately after the integration of the 
completed flight instrument onto the spacecraft platform. In the 
following sections, the instrument visible channels’ band 
responsivities from testing at the instrument manufacturer’s facility 


(using the conventional two-step approach) are compared with the 
Special Test at the spacecraft manufacturer’s facility (using the laser- 
based approach). Detector-to-detector differences, optical cross-talk 
and the uncertainties in the determination of the VIIRS sensor visible 
and near infrared band responsivities are presented. The comparison 
of these two sets of band responsivities engenders a discussion of the 
OOB components of the VIIRS spectral responses and their effects on 
the dependence of the outputs of the bands to the spectral 
distributions of the radiances viewed. 


2. Traveling SIRCUS at BALL Aerospace 


The SIRCUS facility at NIST has been described in detail 
previously [6] with the travelling or T-SIRCUS being used during VIIRS 
sensor testing at BATC. The T-SIRCUS setup provided monochromatic 
radiance from a 76.2 cm diameter SIS equipped with a 25.4 cm 
diameter aperture. A schematic diagram of the NIST T-SIRCUS is 
shown schematically in Fig. 3. Two calibrated Gershun-tube radiance 
meters were used at BATC to provide a calibration of two silicon 
photodiodes mounted on the SIS wall. The calibration related the 
monitor signal from the photodiodes to the sphere radiance. The 
sphere-mounted photodiodes were calibrated pre- and post-VIIRS 
measurements using the NIST working standard radiance meters in 
the BATC high bay outside the clean room. During VIIRS 
measurements of the sphere, the radiance was determined solely by 
the monitor detectors. Before and after the measurement program at 
BATC, the Gershun-tube radiance meters were calibrated at NIST. 
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Figure 3. Schematic diagram of NIST’s T-SIRCUS. 


The tunable laser system remained in the high bay area outside the 
clean room. The integrating sphere was aligned to the VIIRS sensor 
nadir (Earth-view) port in the BATC clean room (Fig. 4). Radiant flux 
from the laser sources was coupled to the SIS using steel-jacketed 
200 um core diameter silica-silica optical fiber. A beamsplitter in the 
optical path sent a small portion of the laser radiation into a 
wavemeter that measured the wavelength of the radiation. A 
remotely controlled electronic shutter in the optical path enabled 
ambient signal levels to be routinely acquired. Room lights were kept 
off during testing. 

During testing, the VIIRS telescope was stationary (non-rotating) 
and VIIRS continuously acquired data at a fixed integration time. T- 
SIRCUS data sets, consisting ofa timestamp, the laser wavelength, and 
the monitor signals, were obtained concurrently with VIIRS 
acquisitions. The shutter was closed while the wavelength was 
changed, providing a beginning and end of wavelength ambient signal 
in the data streams. These ambient signals combined with timestamps 
in the data sets facilitated the proper merging of T-SIRCUS and VIIRS 
data sets. Because the VIIRS integration time was constant, the 
instrument output is recorded as a digital number (DN) rather than a 
DN/unit time. The DNs from the dark periods were subtracted to 
provide the net DNs for each radiance level from the sphere. 


Equation 1 is the measurement equation. S(DN) is the VIIRS signal 
in digital numbers, DN, while viewing the laser-illuminated T-SIRCUS 
integrating sphere, R(A) is the VIIRS responsivity at the measurement 
wavelength A, L(A) is the T-SIRCUS integrating sphere radiance at the 
measurement wavelength, and 6(A - Ao) is the laser lineshape centered 
at the peak wavelength, Ao, represented by a delta function. The VIIRS 





Figure 4. Photograph of the T-SIRCUS SIS in front of the SNPP VIIRS 
sensor nadir doors. 


S(DN) =| R(A)L(A)5(A-A,)d(A) 
S(DN) = R(A,)L(A, )and 


R(A,)=S(DN)/L(A, ) (1) 
net DNs are combined to give an ASR, in DN/(W m~ sr‘), at the laser 
wavelength. A band’s RSR is the peak-normalized ASR. 


3. Historical Approach to Responsivity Calculations 


The instrument manufacturer’s calibration technique for VIIRS 
follows an approach that has been in use for several decades [5, 7]. 
First, the relative spectral responses of the instrument bands are 
determined. Then the instrument views a source with a known 
radiance spectrum such as a lamp-illuminated integrating sphere, and 
the instrument readings for that view are recorded. Finally, for each 
band, the spectral response and the sphere radiance are combined to 
form a band-averaged spectral radiance, which is combined with the 
instrument reading for the band to provide a calibration coefficient. 
Additional quantities of interest include the band-center wavelength 
and bandwidth. 


3.1 Relative Spectral Response Determination 


For the laboratory calibration by the instrument manufacturer, 
Raytheon, the spectral responses for the VIIRS visible and near 
infrared bands (M1 through M7) were measured using the SoMA in 
ambient prior to TVAC testing [8]. The SoMA consists of two MS257 
quarter-meter monochromators used in the subtractive mode, a 
tungsten source or ceramic source (depending on wavelength), and 
collimating optics. The exit beam is collimated by a 30 cm off-axis 
spheroid with an effective focal length (EFL) of 125 cm. The exit slit is 
at tangential focus to give the sharpest image in the scan direction. 
This piece of equipment requires both spectral calibration and relative 
spectral output (RSO) characterization to minimize measurement 
uncertainties. The SoMA produces quasi-monochromatic light with 
bandpass settings of 1 nm or 2.5 nm depending on the wavelengths of 
the bands being measured. For measurements of the radiant flux 
from the SpMA by VIIRS, there is a set of transition optics to focus the 
output of the SpoMA onto the instrument's detectors. The focused 
output of the slit, Fig. 2, covers one band’s detectors. The SoMA needs 
to be re-aligned for each band’s measurements. 

In the Vis and NIR, the source of light for the SpMA is an 
incandescent (tungsten) lamp. It produces a peak radiance near 


1000 nm with decreasing radiance on either side of the peak. This 
causes a spectrally dependent shape in the output from the SpMA’s 
monochromator. To correct for this effect, a reference detector with 
known quantum efficiency also views the flux from the SpMA. This is 
done before and after the VIIRS spectral response measurements. The 
reference detector does not share the transition optics used in the 
measurements by VIIRS, so the detector monitors the relative output 
of the SoMA with wavelength but not the absolute value of that output. 
The SpMA-based spectral responses for each VIIRS band, after 
correction for the spectral shape of the flux from the source, is 
normalized to the peak value for that band’s response to give an RSR. 

A representative SoMA-based RSR spectrum (for VIIRS band 7, 
detector 8) is shown in Fig. 5a. The RSRs are the individual points in 
the figure. The figure includes the wavelength region between the 
band’s 1 % responses relative to the maximum RSR, and the band is 
the compilation of the RSRs. Historically, this wavelength region is 
defined to be the In-Band (IB) region, incorporating the bulk of the 
total integrated band response. Other definitions of the IB wavelength 
region are possible, based for example on the area under the curve 
formed by the individual RSRs in which the IB area incorporates a 
large, fixed percentage of the total area. Each definition is used to 
separate the IB response from the OOB response; the response from 
the OOB component of the total instrument response is typically 
treated as a contamination to the response in the IB response region. 
The Total-Band (TB) response is the summation of the IB response 
and the OOB response. 

A number of quantities are typically used to describe the 
radiometric properties of a particular band. The basic properties ofa 
filter band are often described in terms of its center wavelength and 
spectral width, often called the bandwidth. Combined, these terms 
provide an estimate of the functional wavelength region for the band. 
For quantities such as spectral radiance, spectral irradiance, or the 
band’s center wavelength, the technique of band averaging is 
commonly employed. Band-averaging is defined in Section 3.3 while 
Section 3.2 defines the bandwidth. Both the center wavelength and 
bandwidth can be calculated using the RSRs and their associated 
wavelengths. 
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Figure 5. Spectral Response Measurements for VIIRS band M7, detector 8. 
Measurements are shown over the In-Band response region, where the 
responses are 1 % of full scale and greater: (a) SoMA-based RSR and (b) 
T-SIRCUS-based ASR. 


3.2 Bandwidth 
As the area under the curve, the SpMA-based bandwidth BW, ,, 


can be calculated via trapezoidal integration as 


BWe va” ee a(n -A(n-1)| 
: SE [ase(n) + RSR(n-1)| (2) 


where RSR(n) is the relative spectral response for measurement n, n 
max are the number of measurements, and A is the wavelength (in 
nm) for measurement n. For the SoMA-based RSR measurements, 


individual values in the IB region are separated by a constant 
wavelength ( A ) and the slit widths of the monochromator are 
adjusted to fill that wavelength region at the VIIRS focal plane. 
Because the wavelength spacing AJ is constant, it can be pulled out of 
the summation. Nominal center wavelengths and bandwidths for the 
VIIRS bands studied in this work are given in Table 1; SpMA-derived 
bandwidths are given in Table 2. 


Table 2. Bandwidths and band-center wavelengths for the IB wavelength 
ranges of the moderate resolution VIIRS bands M1 through M7. The values for 
the laboratory measurements come from the SpMA. The 4" and 7® columns 
show the differences, in nm, from T-SIRCUS. 
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Bandwidth SpMA SIRCUS Wavelength 
SpMA SIRCUS Difference Band-Center Band-Center Difference 

VIUIRS Bandwidth | Bandwidth | SpMA-SIRCUS Wavelength Wavelength SpMA-SIRCUS 
Band (nm) (nm) (nm) (nm) (nm) (nm) 
M1 18.53 19.15 -0.62 410.84 410.75 0.09 
M2 14.02 1408 |  -0.06 =| 443.63 443.69 -0.06 
M3 18.87 18.80 0.07 486.13 486.34 -0.21 
M4 19.82 20.04 -0.22 550.94 550.78 0.16 
M5 19.61 19.44 0.17 671.49 671.59 -0.10 











3.3 Band-Averaged Center Wavelength and Spectral Radiance 
Band-averaging uses the spectral shape of the band to create a 
weighted mean ofa particular quantity. For a set of measurements, 

the band average of quantity X ,, is calculated as 





i X (n-1)RSR(n 2 ]A)-20 1] 


NX — n=2 2 


n=2 





, (3) 


where Xaa is the quantity after averaging, and X(n) is the quantity for 
measurement n at wavelength 4 (7) . In this application, X, 


represents either wavelength or sphere radiance. Any constant 
multipliers for these terms can be moved outside the summations and 
cancelled out in the ratio if they are in both the numerator and the 
denominator. The RSR has been measured with a constant 
wavelength spacing. The calibration source radiance is known and is 
often interpolated to the same wavelength spacing. Eq. 3 can then be 
expressed as 


nmax 


; d, | X (n)RSR(n)+ X (n-1)RSR(n-1) | 


X pa nmax 
>. | RSR(n)+ RSR(n-1) | 
a . (4) 
For the SoMA-based band-averaged center wavelength, 
Ana (SpMA) , A(n) is substituted for X(n) in Eq. 4, 
S*[a(n) RSR(n) + A(n-1)RSR(n-1) | 
Ay, ( SpMA) = 2 —___________________ (5) 
>) | RSR(n) + RSR(n-1) | 
n=2 


and for the calibration source band-averaged spectral radiance L;" 
the spectral radiance of the source, 1°” (n) , is substituted for X(n) in 
Eq. 4, 


M6 14.45 14.36 | 0.09 745.41 745.59 -0.18 
862.00 862.10 -0.10 


S"[ 1S (n) RSR(n) + 1°" (n =1) RSR(n 1) | 
L" (SpMA) = =2@—$_£27,—H?____________ 


nmax 


2 | RSR(n) + RSR(n—-1)| 
(6) 


The summation extends over the full wavelength range where 
there is measurable response from the band. Note that these band- 
averaged quantities are TB quantities, as opposed to IB quantities, 
where the summation is restricted to a defined spectral range 
corresponding to the +1 % relative response wavelengths. 


3.4 Band Responsivity Calculations 


For traditional responsivity calculations, the spectral radiance 
values at individual wavelengths comes from the radiance spectrum of 
a calibrated lamp-illuminated integrating sphere and its weighted 
mean is the band-averaged spectral radiance. Because the summation 
is over the full spectral response region, the resulting quantity is the 
TB Band-Averaged Responsivity. The band-averaged responsivity 

R,, is defined to be the instrument output DN acquired while 


measuring the calibration source divided by the calculated band- 


L™" (SpMA 
averaged spectral radiance ™ (Sp ) 


DN 


R,,, (SpMA) = —————. 
ns (SPMA) Lys (SpMA) 


(7) 
Since the radiance source for the instrument manufacturer's 
calibration, the SIS100, is calibrated in VIIRS units, (W m~ sr? pn), 
R,, has units of DN/(W m* sr? um"). 
For this approach, there are two principal uncertainty components 
in the calibration of the VIIRS bands: the uncertainty in the RSRs from 


the SoMA measurements and the uncertainty in the spectral radiance 
from the calibration source, in this case the SIS100. 


4. TUNABLE LASER (T-SIRCUS) APPROACH TO 
RESPONSIVITY CALCULATIONS 


The laser-based measurements give an ASR at each measurement 
wavelength. A representative ASR spectrum (for VIIRS band M7, 
detector 8) is shown in Fig. 5b. Measurements over the IB response 
are shown, corresponding to the spectral range shown in Fig. 5a for 
the SoMA-based RSR measurements. The laser sources used in the 
calibration of the VIIRS bands are essentially delta functions and the 
spacing of the wavelengths in the ASR measurements can be 
arbitrarily small, as opposed to the SoMA approach where the step 
size is matched to the monochromator bandpass. As shown in Fig. 5b, 
the spacing of the laser-based ASR measurements near the central 
peak is sufficiently small to show any structure within the band 
response and to allow for high accuracy band integration. For direct 
comparison with the SpMA-based results, an ASR can be converted to 
a normalized RSR via division by the maximum ASR, a: 


1 
RSRyiecus (n) = oy to Rsiacus (1) , (8) 


where « has units of DN/(W m-2 sr-1). For the ASR spectrum in 

Fig. 5b, for example, the maximum ASR occurs at 850.3 nm and a hasa 
value of 2860.5 DN/(W m-2 sr-1). In the discussion in Sections 4.1 to 
4.3, itis understood that ASRs and RSRs refer to T-SIRCUS-based 
measurements. 


4.1 Laser-based Bandwidth Derivation 


A combination of Eqs. 2 and 8 provides the T-SIRCUS-based 
bandwidth BWo cus : 


BWoyoss = 5 |e eee es JE) <a) 
e (9) 


nmax 


= = ¥[(45K(0) + ASR(n-1))], 


a n=2 


where (1/a) is applied to each ASR in the equation, maintaining the 
units from Eq. 2. 


4.2 Calculated Band-Integrated Responsivity and Band-Averaged 
Center Wavelength 


With the laser-based approach, there is no broad-band calibration 
source and the band-integrated responsivity R,, is often calculated in 


place of the band-averaged responsivity R,,. R,, can be obtained 
through trapezoidal summation, 


p= S| ROAD a) ao], (10) 


n=2 


where ASR(n), the ordinate, is the monochromatic ASR for 
measurement n, and A(n), the abscissa, is the wavelength for that 
measurement. 

The combination of Eqs. 4 and 8 provides the T-SIRCUS-based 
band-averaged center wavelength /,,(SIRCUS) , 


nmax 


> [A(n) ASR(n) + 2(n-1) ASR(n-1)] 


A,, (SIRCUS ) = = (11) 


nmax 


> [ASR (1) + ASR(n-1)] 

Eq’s.9 and 11 allowa direct comparison of bandwidths and band- 
center wavelengths derived from T-SIRCUS-based laboratory 
measurements with those quantities derived from SpMA 
measurements. 


4.3 Spectral Radiance from a Known Source 

For the traditional instrument calibration, the known spectral 
radiance spectrum from a lamp-illuminated integrating sphere L™ is 
combined with the VIIRS instrument band’s ASR to provide a 
calibration band-averaged spectral radiance, Le .A combination of 


Eqns. 4 and 8 allows the calculation of L;"" , again assuming constant 
wavelength spacing: 


Se (pane) (Naso) 


L,” (SIRCUS ) = = (12) 


nmax 


> [ASR (n) + ASR(n -1)] 


n=2 


5. CONVERSION FROM TOTAL-BAND TO IN-BAND 
VALUES 


The band-averaged quantities derived in Sections 3 and 4 are all TB 
quantities because the summations run over the full measurement 
range, which should coincide with the spectral range of finite response 
of the instrument band. The desired values from instrument bands 
are most often those near the central peak of each band’s response, or 
from the IB region. 

For IB quantities (band-averaged center wavelength, band- 
averaged radiance, band-averaged responsivity), the summation is 
over a finite spectral region. For VIIRS, the IB region is defined to be 
the spectral region between the +1 % relative responsivity values. A 
measured response in the OOB region contributes a spectrally 
dependent error to IB quantities. To eliminate OOB biases in TB 
quantities, often the TB quantities are converted to give their IB 


counterparts by subtracting the OOB contribution from the TB 
quantity. 

For the T-SIRCUS-based band-integrated responsivities, which are 
calculated as the sum of individual ASRs, the process is the same for 
calculating R,, as for R;; with reduced limits on the summation. The 


SpMA-based calculation of the Total-Band band-averaged 
responsivity R;; (SpMA) in Eq. 7 includes the laboratory calculation 


of the instrument output in DNs. For the SoMA-based conversion from 
a TB value to an IB value, the solution involves the derivation of the 
RSR-to-ASR scaling factor a based on measurements of a source with 
known spectral radiance, in this case the SIS100. The measurement 
equations for a VIIRS band-measurement of a calibration source are: 


2 


ae 5 LE" (n) ASR(n) + L™ (n= 1) ASR(n- Dat) oe 
™ (13) 


- ale (n) RSR(n) + L™ (n-1) RSR(n-1) |[A(n) + 4(n-1)] 


and 


n max 


a =2DN / Se (n) RSR(n) + L (n=1) RSR(n - D|LG@) +A(n- 1] 


(14) 
Once ais known, R;, (SpMA) can be calculated the same way that 


R;, (SIRCUS ) is calculated. 


6. COMPARISONS 


The basic values in the instrument characterization, bandwidth and 
band-averaged center wavelength are described in Sections 3 along 
with the band-averaged responsivity, the basic quantity in the 
instrument calibration. A comparison of instrument characterization 
results using the two methodologies is presented here. RSR 
comparisons are made at individual wavelengths; for these 
comparisons, the T-SIRCUS ASRs are converted to RSRs using Eq. 8 
and then linearly interpolated to the SoMA RSR wavelengths. 


6.1 Comparison of IB RSRs 


Comparisons of the IB RSR between the SpMA-based and the T- 
SIRCUS-based derivations are shown in Figs. 6-12. In these figures, 
the SpMA and interpolated T-SIRCUS RSRs are shown in the a) panels 
while the differences (SpMA from T-SIRCUS) are shown in the b) 
panels. The b) panels also provide an outline of the T-SIRCUS RSRs as 
a visual reference for the locations of the band edges. The T-SIRCUS 
ASRs have been normalized to unity and interpolated to the SoMA 
wavelengths. Since T-SIRCUS provides a greater number of 
measurements per unit wavelength, the interpolation has removed 
details from the set of T-SIRCUS RSRs. The features in the response 
differences at the band edges exhibiting a rapid change in in response 
in these panels could arise from wavelength mismatches between the 
measurements by T-SIRCUS and by the SpMA,; the differences could 
also be a reflection of the finite 1 nm or 2.5 nm bandpasses set in the 
SpMA monochromator. There is a 0.5 nm uncertainty in the 
wavelength characterization of the VIIRS bands using the SpMA, 
determined from views of atomic emission line spectra. For the 
wavelengths in T-SIRCUS, the uncertainties are less than 0.01 nm for 
the VIIRS measurements, except for wavelengths near 500 nm where 
the uncertainty is 0.2 nm and the laser power is low. 

For band M1 in Fig. 5a, the T-SIRCUS-measured response has a 
peak near 413 nm, and the SpMA peak is near 417 nm. This may 
result from a downward drift in the output of the SoMA lamp in the 
blue during the interval between the lamp’s characterization (as 
described in Section 3.1) and the measurement of the band’s RSR. 
Instabilities in the SoMA lamps were problems for the laboratory 
characterization of VIIRS. During some RSR measurements, the SoMA 
lamps experienced failures, and some failed before their second 


relative spectral output (RSO) characterization following the RSR 
measurements. Although it is not known with certainty, output 
instability might apply to the lamp used for band M1. In general, the 
effects of SoMA lamp instabilities remain difficult to quantify. 


6.2 Bandwidth, Center Wavelength, and_ Responsivity 
Comparisons 


The IB bandwidths and band-center wavelengths for the SoMA and 
T-SIRCUS data sets are listed in Table 2 and the band responsivities in 
Table 3. The summations for the T-SIRCUS calculations for these 
comparisons use the full set of laboratory measurements and are not 
interpolated to the SoMA wavelengths as in the individual wavelength 
comparisons above. For bandwidth and band-averaged center 
wavelength, the differences from T-SIRCUS (in nm) are also listed in 
Table 2; and for responsivity, the percent differences from T-SIRCUS 
are listed in Table 3 and are displayed in Fig. 13. For the center 
wavelengths, the average differences are well within the 0.5 nm 
wavelength characterization uncertainty for the SpMA. For the IB 
responsivities, the differences (SoMA from T-SIRCUS) average 2.2 % 
(+1.7 %, k=1). This isa measure of the agreement in the calibrations 
from the two techniques. 


6.3 Comparisons for the TB RSRs 


For the seven VIIRS bands, the SpMA and T-SIRCUS RSRs are 
shown over their entire sets of measurement wavelengths in Figs. 14- 
17. The ordinates are logarithmic, covering six orders of magnitude to 
improve visualization of the OOB response. For T-SIRCUS, the RSRs 
are calculated from the ASR measurements via division bya The full 
data set is shown with individual data points at their measured 
wavelengths. Gaps in the T-SIRCUS measurements near 550 nm (Fig. 
14 in particular) stem from low output power from the T-SIRCUS laser 
in this wavelength range. 

The TB values for these figures are calculated in the same manner 
as the IB values in Section 6.1. The bandwidths and center 
wavelengths for the two data sets are listed in Table 4, and the 
responsivities are listed in Table 5. The sets of differences are 
displayed in Fig. 18. There is considerably more scatter in the TB 
wavelength differences (Fig. 18a) relative to the IB differences (Fig. 
13a). This is due in part to the limited widths of the In-Band responses 
relative to the overall wavelength range of the measurements (Fig 14, 
for example). Asa result, the differences in the OOB regions become 
significant. In addition, the center wavelengths are calculated using 
weighted means, with the wavelengths as the weighting factors (see 
Eq. 3). For bands M1 and M2, in particular, the wavelengths in the red 
and near-infrared are a factor of two greater than those for the IB 
responses, amplifying the contribution of the OOB differences. For the 
bandwidths, on the other hand, the TB differences are in much closer 
agreement with their IB counterparts. For the TB responsivities, the 
differences (SpMA from T-SIRCUS) average 2.0 % (+1.6 %,, k=1), in 
close agreement with the IB differences. 

The maximum band responses (a’s) are listed in Table 6. For T- 
SIRCUS, these values are determined from each set of measurements 
by inspection. For the manufacturer’s measurements, the maximum 
ASRs come from the division of the responsivities by the bandwidths 
via Eq. 14. For the maximum ASRs, the differences (SpMA from T- 
SIRCUS) average 2.6 % (+2.0 %, k=1), in agreement with both the IB 
and the TB responsivity differences. 


6.4. Conversions from Total-Band to In-Band Responses. 


The desired values from the instrument are those near the central 
peak of each band’s response, that is, within the IB wavelength region. 
To obtain them, itis necessary to convert the TB values into their IB 
counterparts. For bandwidth and center wavelength responsivity, this 
can be done with the difference (in nm) of the IB values from their TB 
counterparts. The corrections for bandwidth are listed in Table 7 and 


those for center wavelength in Table 8. The corrections are plotted in 
the a) and b) panels of Fig. 19, respectively. 

The nominal center wavelength for Band M1 is 410 nm (see 
Table 1). The IB values are within 1 nm of nominal (see Table 2). 
However, the TB center wavelengths have values near 421 nm (see 
Table 3). These differences come from the relatively large OOB 
response for the band (Fig. 14a), plus the sensitivities in the 
wavelength calculation described in Section 6.2. In a similar manner, 
there are corrections close to 2 nm for bands M2 and M3. The relative 
differences in the SoMA and T-SIRCUS wavelength corrections for 
band M2 in Fig 19b can also be seen in the Total-Band wavelength 
differences in Fig. 18a. They result from the greater OOB response 
from T-SIRCUS in Fig. 14b. A possible explanation for this difference is 
discussed in the Crosstalk section below.The bandwidth corrections 
(Fig. 18a) are less than 1 nm in all cases. For band M4, the bandwidth 


Table 3. Band responsivities for the IB wavelength ranges of VIIRS 
bands M1 through M7. The responsivities for the laboratory 
measurements come from the SpMA and the SIS100. The comparisons 
show the differences between the SpMA and T-SIRCUS. 


Responsivity Difference 
VIIRS SpMA Band SIRCUS Band SpMA-SIRCUS 
ia Responsivity' Responsivity' (%) 


‘Units: DN/(W m? sr! um’) 





Table 4. Bandwidths and center wavelengths for the TB wavelength 
ranges of VIIRS bands M1 through M7. The responsivities for the 

laboratory measurements come from the SpMA and the SIS100. The 

comparisons show the differences between the SpMA and T-SIRCUS. 


SpMA SIRCUS | Bandwidth SpMA SIRCUS Wavelength 
Bandwidth | Bandwidth | Difference | Wavelength | Wavelength Difference 
(nm) (nm) (nm) (nm) (nm) (nm) 

19.08 19.70 -0.62 421.49 420.65 0.84 
-0.18 443.92 445.75 1.83 
19.06 19.08 -0.02 488.74 489.46 -0.72 
M4 20.67 20.89 -0.22 551.58 552.07 -0.49 
M5 20.18 20.00 0.18 670.24 671.01 -0.77 
14.67 14.59 0.08 744.62 745.10 -0.48 


38.89 38.68 0.21 861.82 861.74 0.08 





























Table 5. Responsivities for the TB wavelength ranges of VIIRS bands M1 
through M7. The responsivities for the laboratory measurements come from 
the SpMA and the SIS100. The differences between the SoMA and T-SIRCUS- 
based calibration are given in column 4. 


VIIRS SpMA/SIS SIRCUS Responsivity 
Band Responsivity’ Responsivity’ Difference (%) 


M3 27.568 21,053 
M4 38.427 36.888 





M5 52.283 51.129 
M6 86.195 85.410 








M7 111.892 110.639 
‘Units: DN/(W m” sr! um’) 











correction has the greatest magnitude. For this band, the RSRs have 
values greater than 0.003 for a wavelength range of about 100 nm 
(Fig. 15b). It has the greatest accumulated OOB response for the VIIRS 
bands. 


Table 6. Maximum Absolute Spectral Responses (ASRs) for the Total-Band 
wavelength ranges of VIIRS bands M1 through M7 using the SpMA/SIS and T- 
SIRCUS approaches. 

SpMA/SIS SIRCUS Difference 
Maximum Maximum (%) 
ASR! ASR 


1046.5 1014.7 
1741.8 1656.2 





1 -2 _-1 
sr 


Units: DN/(Wm ) 


For responsivity, the correction factor (Table 9) is multiplicative, 
that is, the In-Band value divided by the Total-Band value. This 
conversion factor is dimensionless. The correction factors for the 
bands are shown in Fig. 19c. The pattern in this panel is similar to the 
pattern for bandwidth in panel a) of the same figure. The RSRs and 
the ASRs for each band are related by the maximum ASR (a, see Eq. 8), 
which is constant for the band. Asa result, for each band the 
conversion factors for bandwidth and responsivity are identical in 
multiplicative terms. In other words, the multiplicative corrections in 
Table 9 can be applied to the bandwidths in Table 7 as well. 

The difference between the TB and IB responsivity responses 
comes from the OOB response. For example, the conversion factor for 
the band M1 SpMA responsivities in Table 9 is 0.9713. For this band, 
the IB portion is 97.13 % of the TB and the OOB is the other 2.87 %. 
For the band M1 T-SIRCUS responsivities, the partitioning of the Total- 
Band response is 97.17 % In-Band and 2.83 % OOB (see Table 9). For 
band M1, the OOB contributions for the SoMA and for T-SIRCUS differ 
by 0.04 %, a practically insignificant amount. For bands M4, M5, and 
M6, the differences between the SpMA and T-SIRCUS OOB 
contributions to the TB are equally small, 0.05 %, 0.02 % and 0.06 %, 
respectively (see Table 9). In other words, the OOB differences 
between their SpMA and T-SIRCUS responses are practically 
insignificant for these bands as well. Thus, for band M1, the SoMA/T- 
SIRCUS differences from 450 nm to 600 nm that appear substantial in 
the logarithmic plot have almost no effect on the conversion from the 
TB to the IB. These differences have been attributed to very weak or 
to missing monochromatic radiances caused by problems with the 
lasers used for the T-SIRCUS measurements. For bands M4, M5, and 
M6, the noise floor in the T-SIRCUS measurements for wavelengths 
below about 440 nm (see Figs. 15b, 16a, and 16b) is nota significant 
contributor to the TB/IB conversion either. 

For bands M2, M3, and, M7, the SoMA and T-SIRCUS conversion 
factors for responsivity differ by more significant amounts (see Table 
9 and Fig. 19c). For each of these conversion factors, the SpMA value is 
closer to unity than its T-SIRCUS counterpart, indicating that the SoMA 
OOB contributions are smaller (Table 9). For band M2, the difference 
is 0.85 %, and for bands M3 and M7, the differences are 0.43 % and 
0.17 %, respectively. As explained in the following section, these 
differences are attributed to crosstalk between bands in the 
instrument, an effect that the SoMA under-represents. 


7. CROSSTALK 


Unexpected features were observed in the VIIRS sensor response 
during spatial response testing of the Engineering Design Unit (EDU). 
Examination of the phenomenon led to its identification as crosstalk, 
or communication between bands. Multiple types of crosstalk were 
identified, each with its own pathway for band-to-band 
communication. The two main types crosstalk were called dynamic 
and static crosstalk. 

Dynamic crosstalk is an electronic effect where a change in the 
photo-current in one band invokes a response in another 
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Figure 6. Comparison of RSRs from the SpMA and T-SIRCUS for VIIRS band M1, detector 8. This comparison shows differences (SpMA from T- 
SIRCUS) at the SoMA wavelengths. a) T-SIRCUSand SpMARSRs; b) Response Differences (SpMA from T-SIRCUS). 
An outline of the T-SIRCUS RSR is provided as a visual reference for the locations of the band edges. 
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Figure 7. Comparison of RSRs from the SpMA and T-SIRCUS for VIIRS band M2, detector 8. This comparison shows differences (SpMA from T- 
SIRCUS) at the SoMA wavelengths. a) T-SIRCUS and SpMA RSRs; b) Response Differences (SpMA from T-SIRCUS). 
An outline of the T-SIRCUS RSR is provided as a visual reference for the locations of the band edges. 
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Figure 8. Comparison of RSRs from the SpMA and T-SIRCUS for VIIRS band M3, detector 8. This comparison shows differences (SoMA 
from T-SIRCUS) at the SoMA wavelengths. a) T-SIRCUS and SpMA RSRs. b) Response Differences (SpMA from T-SIRCUS). 
An outline of the T-SIRCUS RSR is provided as a visual reference for the locations of the band edges. 
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Figure 9. Comparison of RSRs from the SpMA and T-SIRCUS for VIIRS band M4, detector 8. This comparison shows differences (SpMA from T- 
SIRCUS) at the SoMA wavelengths. a) T-SIRCUS and SpMA RSRs. b) Response Differences (SpMA from T-SIRCUS). 
An outline of the T-SIRCUS RSR is provided as a visual reference for the locations of the band edges. 
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Figure 10. Comparison of RSRs from _ SpMA ee T -SIRCUS vr VIIRS band M5, detector 8. This comparison shows differences (SoMA 
from T-SIRCUS) at the SoMA wavelengths. a) T-SIRCUS and SpMA RSRs. b) Response Differences (SpMA from T-SIRCUS). 
An outline of the T-SIRCUS RSR is provided as a visual reference for the locations of the band edges. 
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Figure 11. Comparison of RSRs from the SpMA and T-SIRCUS for VIIRS band M6, detector 8. This comparison shows differences (SoMA 


from T-SIRCUS) at the SoMA wavelengths. a) T-SIRCUS and SpMA RSRs. b) Response Differences (SpMA from T-SIRCUS). 
An outline of the T-SIRCUS RSR is provided as a visual reference for the locations of the band edges. 
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Figure 12. Comparison of RSRs from the SpMA and T-SIRCUS for VIIRS band M7, detector 8. This comparison shows differences (SoMA 
from T-SIRCUS) at the SpoMA wavelengths. a) T-SIRCUS and SpMA RSRs. b) Response Differences (SpMA from T-SIRCUS). 
An outline of the T-SIRCUS RSR is provided as a visual reference for the locations of the band edges. 
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Figure 13. Differences between T-SIRCUS and the SpMA-based derivations of the In-Band bandwidths, center wavelengths, and responsivities of 
the laboratory-based calibration by the instrument manufacturer. a) Differences of the SoMA-based bandwidths and band center wavelengths 
from T-SIRCUS. b) Differences of the SoMA-based responsivities from T-SIRCUS. 
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Figure 14. Total-Band RSRs for bands M1 and M2. The open squares come from SpMA measurements; the closed circles from 
T-SIRCUS. 
a) RSRs for VIIRS band M1; b) RSRs for VIIRS band M2. 
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Figure 15. Total-Band RSRs for bands M3 and M4. The open squares come from SpMA measurements; the closed circles from T-SIRCUS. 
a) RSRs for VIIRS band M3,; b) RSRs for VIIRS band M4. 
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Figure 16. Total-Band RSRs for bands M5 and M6. The open squares come from SpMA measurements; the closed circles from T-SIRCUS. 
a) Total-band ASRs for VIIRS band M5 &b) RSRs for VIIRS band M6. 
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Figure 17. Total-Band RSRs for band M7. The open squares come 
from SpMA measurements; the closed circles from T-SIRCUS. 
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Figure 18. Differences from T-SIRCUS of the Total-Band bandwidths, 
center wavelengths, and responsivities of the laboratory-based 
calibration by the instrument manufacturer. 

a) Differences of the SoMA-based bandwidths, and center wavelengths 
from T-SIRCUS (in nm). b) Differences of the SoMA/SIS100-based 
responsivities from T-SIRCUS (in percent). 


Figure 19. Correction factors changing Total-Band bandwidths and 
responsivities to their In-Band counterparts. 
a) Correction factors for bandwidth. 
b) Correction factors for band-center center wavelength. 
c) Correction factors for responsivity. 


band. Although this was initially the strongest effect observed, it was 
the easiest to mitigate by application of bonding wires to reduce a 
higher than expected sheet resistance in the common bias electrode of 
the detector array. Subsequently, dynamic crosstalk became small 
relative to static crosstalk. 

Two types of static crosstalk were identified from the VisNIR focal 
plane: electrical and optical. Electrical crosstalk generally occurred 
systematically from point to point across the detector array. This 
phenomenon was also significant in the EDU and led to re-working of 
the electronics to minimize the effect. Static electrical crosstalk was 
found to be small in the flight instrument relative to other static 
effects. Optical static crosstalk, or simply optical crosstalk, occurs 
between a source of the scattered light (sending band) and a recipient 
(receiving band). In the VIIRS EDU, optical crosstalk stemmed from 
incident light in one band scattered at the surface or within the 
Integrated Filter Array (IFA) and then falling into a different band. It 
was noted that, in most cases, the sending sources were predominant 
over many narrow windows in the NIR region, regardless of the 
center wavelength of the sender band. Bench test measurements 
from the SNPP VIIRS flight IFA suggested that scattering from the filter 
array to be greater than anticipated [9], leading to concerns about the 
quality of ocean color data products [10] and other data products that 
are sensitive to radiometric accuracy. However, the more significant 
effects of optical crosstalk were limited to adjacent bands on the focal 
plane, and may have been artifacts of the bench experiment 
employed. Meanwhile, the IFA was returned the focal plane filter 
bezel with the opposite side facing incoming light. Based on bench 
tests of the EDU IFA, it was decided that flipping the IFA would reduce 
the spread of light between the IFA and the detector array 
beneath. This was effective, but with a trade with significantly 
increased out-of-band response. 

T-SIRCUS was applied at the sensor level with the final IFA 
orientation. The T-SIRCUS input radiant flux floods the entire focal 
plane in the same manner as measurements of upwelling light from 
the Earth on-orbit and all forms of channel-to-channel crosstalk is 
automatically included in the measurements, along with out-of-band 
response. For VIIRS bands M1 and M3, the T-SIRCUS and SPMA- 
derived RSRs show series of response maxima at wavelengths greater 
than about 600 nm (see figures 14a and 15a). These effects were 
attributed to optical crosstalk. For VIIRS bands M2, M5, and M7, the T- 
SIRCUS-derived RSRs show responses from wavelengths greater than 
about 600 nm not found in the SoMA-derived RSRs. These differences 
were attributed to static electrical crosstalk, which are distinguished 
by their clear outline of the in-band pass of the sending band. Figures 
16a and 16b show a similar crosstalk effect near 860 nm for bands M5 
and M6, respectively, while Fig. 17 shows band M7 crosstalk effects 
from bands M5 and M6 near 680 nm and 745 nm, 
respectively. Crosstalk features are explicitly labeled in the figure. For 
band M2 (Fig. 14b), there are several crosstalk contributions in the T- 
SIRCUS RSRs at wavelengths in the red and near infrared. 

Simulation of crosstalk effects on SNPP VIIRS ocean color band 
radiometry at the Top-of-Atmosphere (TOA) demonstrated that the 
magnitude of the crosstalk effect ranged up to few tenths ofa 
percent [11]. This was much reduced compared to the earlier tests of 
the EDU and more manageable for radiometrically sensitive 
products. The increased out-of-band, conversely, was up to a few 
percent in the worst case. However, this was still in line with crosstalk 
magnitudes observed in heritage instruments such as MODIS and 
SeaWiFS [12], and considered correctable. The SpMA provided 
important information regarding the spatial characteristics of 
crosstalk. T-SIRCUS, however, provided more realistic information 
regarding the behavior of the combined types of crosstalk and out-of- 
band response for relatively uniform targets in TOA imagery. The 
inclusion of any residual crosstalk response would be included in the 


T-SIRCUS RSR, effecting a correction for low contrast targets in TOA 


imagery. 


Table 7. TB and IB bandwidths from the instrument manufacturer’s 
calibrationand from T-SIRCUS. For both, the conversion factors are 
Total-Band bandwidths minus the In-Band values. 



































VIIRS SpMA SpMA Bandwidth SIRCUS SIRCUS Bandwidth 
Band In-Band Total-Band Conversion In-Band Total-Band Conversion 
Bandwidth Bandwidth Factor Bandwidth Bandwidth Factor 
(nm) (nm) (nm) (nm) (nm) (nm) 
Mi 18.53 19.08 -0.55 19.15 19.70 -0.55 
M2 14.02 14.07 -0.05 14.08 14.25 -0.17 
M3 18.87 19.06 -0.19 18.80 19.08 -0.28 
M4 19.82 20.67 -0.85 20.04 20.89 -0.85 
M5 19.61 20.18 -0.57 19.44 20.00 -0.56 
M6 14.45 14.67 -0.22 14.36 14.59 -0.23 
M7 38.74 38.89 -0.15 38.46 38.68 -0.22 


Table 8. TB and IB center wavelengths from the instrument 
manufacturer’s calibration and from T-SIRCUS. For both, the conversion 
factors are TB center wavelengths minus the In-Band values. 





VIIRS 
Band 


SpMA 
In-Band 
Wavelength 
(nm) 


SpMA 
Total-Band 
Wavelength 
(nm) 


Wavelength 
Conversion 
Factor (nm) 


SIRCUS 
In-Band 
Wavelength 
(nm) 


SIRCUS 
Total-Band 
Wavelength 

(nm) 


Wavelength 
Conversion 
Factor (nm) 





Mi 


410.84 


421.49 


-10.65 


410.75 


420.65 


-9.90 





M2 


443.63 


443.92 


-0.29 


443.69 


445.75 


-2.06 





M3 


486.13 


488.74 


-2.61 


486.34 


489.46 


-3.12 





M4 
M5 


550.94 
671.49 


551.58 
670.24 


-0.64 
1.25 


550.78 
671.59 


552.07 
671.01 


-1.29 
0.58 





M6 


745.41 


744.62 


0.79 


745.59 


745.10 


0.49 








M7 





862.00 





861.82 








0.18 





862.10 





861.74 





0.36 





Table 9. TB and IB responsivities from the instrument manufacturer’ s 
calibration and from T-SIRCUS. For both, the conversion factors are the 
In-Band responsivities divided by the TB values. 























Responsivity Responsivity 

SpMA/SIS SpMA/SIS Conversion SIRCUS SIRCUS Conversion 
VIIRS In-Band Total-Band Factor In-Band Total-Band Factor 

Band Responsivity Responsivity' | (dimensionless) | _Responsivity' Responsivity' | (dimensionless) 

M1 19.392 19.966 0.9713 19.427 19.993 0.9717 
M2 24.421 24.505 0.9966 23.321 23.600 0.9881 
M3 27.286 27.568 0.9898 26.660 27.053 0.9855 
M4 36.855 38.427 0.9591 35.399 36.888 0.9596 
MS 50.815 52.283 0.9719 49.701 51.129 0.9721 
M6 84.892 86.195 0.9849 84.068 85.410 0.9843 
M7 111.470 111.892 0.9962 110.028 110.639 0.9945 





























‘Units: DN/CW m? sr! um’) 


8. DETECTOR-TO-DETECTOR DIFFERENCES 


The calculations in the previous sections have involved a single 
detector, detector 8, of each band in the focal plane. On orbit, each 
scan of the VIIRS mirror sweeps out 16 detectors in the along-track 
direction so there are 15 additional detectors per band in the focal 
plane. Detector 16 in the first scan is adjacent to detector 1 in the 
second scan, and so forth for succeeding scans. Striping occurs in the 
processed imagery if the differences in detector responsivities are not 
accounted for while differences in band-center wavelengths will 
impact data products if not considered. This section gives the 
detector-to-detector band-center wavelength and band-averaged 
responsivity differences for the set of band M1 detectors in the focal 
plane as a representative example of effects seen in all bands. 

Spectral distortion in an imaging spectrograph is called smile. We 
keep the same terminology to refer to detector-to-detector band- 
center wavelength differences in the focal plane. The band-center 
wavelengths are calculated using Eq. 5 for the set of 16 SoMA-based 
RSRs and Eq. 11 for the set of 16 T-SIRCUS-based RSRs. The relative 
band-center wavelength change across the focal plane is derived 
relative to detector 8 by subtracting each detector’s band center 
wavelength from detector 8’s band center wavelength. For both sets 
of RSRs, the results are shown in Fig. 20a. For the SoMA-based band- 
center wavelengths, the wavelength differences increase in a roughly 
linear fashion from about -1.5 nm for detector 1 to almost 3 nm for 
detector 16. The T-SIRCUS-based wavelength differences agree well 
with the SpMA values for detectors 1 to 8. However, for detectors 9 
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to 16, the T-SIRCUS-based band-center wavelength differences flatten, 
with the wavelength for detector 16 about 0.3 nm greater than 
detector 8. 

For responsivity, the individual SoMA-based responsivities are 
calculated using Eq. 7 and the T-SIRCUS-based responsivities via Eq. 
10. The detector-to-detector responsivity differences for band M1 are 
shown in Fig. 20b as the percent difference of each detector from 
detector 8. Similar to the band-center wavelength smile, the SpMA- 
and T-SIRCUS-based differences in Fig. 20b agree well with each other 
from detector 1 to 7. However, from detector 9 to 16 they diverge. 
For detector 16, the SoMA-based responsivity differences from 
detector 8 are 8.2 % while the T-SIRCUS-based differences are 3.4 %. 

If the relative responsivities for the 16 detectors in band M1 are 
correctly determined by the T-SIRCUS-based smile, then the detector 
16 to detector 1 transitions on orbit have a step function in the 
measured radiance equivalent to the difference in responsivity 
between detector 16 and detector 1 of 3.4 %. However, if the SoMA- 
based smile is correct, then each detector 16 to detector 1 transition 
on orbit has a step function in radiance of almost 9 %. This step in 
radiance would result in a significant striping in the on-orbit data set 
in the along-track direction. From the start of SNPP VIIRS on-orbit 
operations, striping has been present for band M1, including a change 
in radiance of 5 % from detector 16 to detector 1. Striping is 
ameliorated on orbit by using the solar-illuminated VIIRS diffuser as a 
source of uniform radiance and normalizing the individual detector 
responsivities to give the same measured radiance from the diffuser. 
Consequently, to date, there has been no use of the T-SIRCUS-based 
detector-to-detector responsivities or spectral smile as part ofa 
striping reduction process. 


9. LABORATORY CALIBRATION UNCERTAINTIES 


The uncertainties in the wavelength characterizations from the two 
techniques were given in Section 6.1. The uncertainties for the 
calibration coefficients (responsivities‘) are provided in this section. 
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Figure 20. Detector-to-detector differences for VIIRS band M1. The 
differences are normalized to detector 8. The SpMA-based differences are 
shown as open squares; the T-SIRCUS-based differences as solid circles. a) 

Center wavelength differences in nm. b) Responsivity differences in 

percent. 


For the historical calibration approach used by the VIIRS 
instrument manufacturer, the calibration coefficients are calculated 
from the RSRs and the SIS100 radiance spectrum (see Eq. 7). For T- 
SIRCUS, the calibration coefficients are calculated from the set of 
individual ASRs (see Eq. 12). 


9.1 Uncertainty in the Historical Calibration Approach 


No known uncertainty budget was developed for the 
manufacturer's radiance calibration of SNPP VIIRS because it was 
designed to operate in reflectance-mode. However, the uncertainties 
of the derived band responsivities are principally a function of the 
uncertainty in the derivation of each band’s RSR and the uncertainty in 
the spectral radiance from the SIS100. To establish the industry-wide 


uncertainty in the radiance from lamp-illuminated SISs, NASA’s Earth 
Observing System (EOS) Project Science Office supported a set of 7 
validation campaigns to satellite sensor calibration facilities and 
measured the radiance of lamp-illuminated SISs at those facilities 
using a set of independently calibrated transfer radiometers [13]. The 
validation campaign ran from 1995 through 2001. From an analysis 
of the resultant data set, the industry-wide mean uncertainty in the 
radiance from lamp-illuminated SISs was + 3 % in the visible to near 
infrared spectral region. 

The lamp-illuminated SIS calibration and disseminated scale 
validation techniques have remained the same over the intervening 
years, so we believe the 3 % standard (k= 1) uncertainty from NASA 
EOSs validation campaign applies to the current uncertainty in the SIS 
radiances and, thus, to the manufacturer’s calibration of SNPP VIIRS. 
This is a minimal uncertainty in the historical SoMA/SIS100 
calibration approach because there are additional sources of 
uncertainty in the approach, including the SpMA-based wavelength 
uncertainty and finite bandpass uncertainty. 


9.2. Uncertainty in the Tunable Laser (T-SIRCUS) Approach 


The uncertainty budget for T-SIRCUS measurements of SNPP VIIRS 
is given in Table 10. The uncertainty budget applies to the In-Band 
portion of the VIIRS responses, not to the Out-of-Band regions where 
the ASRs can be five orders of magnitude or lower below a band’s 
peak response. Individual elements are listed, and their uncertainty 
values for SNPP VIIRS are listed in the third column. The table is 
accurate for the spectral range covered with the exception of the 
region around 550 nm, where the wavelength uncertainty increased 
to 0.2 nm. The final column in Table 10, Future Target, illustrates 
potential T-SIRCUS calibration uncertainties going forward. 

The dominant uncertainty component in the SNPP VIIRS 
measurements is the uncertainty arising from the spatial non- 
uniformity in the T-SIRCUS SIS radiance. It was set to reflect the two 
different calibration coefficients for the monitor photodiodes on the 


Table 10. NIST T-SIRCUS uncertainty budget for SNPP VIIRS. 


[ae [sr 


Reference Detector Stability 
(change! 





SIS wall and was set to an upper limit to be conservative in the 
component's magnitude. The interpolation error arises from the low 
density of points in the calibration of the NIST reference detector, not 
the VIIRS sensor responsivity. The combined standard uncertainty in 
the SNPP VIIRS calibration by the T-SIRCUS is estimated to be 0.28 %; 
the expanded (k = 2) uncertainty is 0.56 %. 

With enhancements to the T-SIRCUS laser system (eliminating the 
lower power and larger wavelength uncertainty in the 500 nm to 
550 nm spectral range), and mapping the SIS spatial radiometric 
uniformity, the potential expanded uncertainties in the T-SIRCUS 
calibration of a VIIRS instrument in the silicon range can be reduced to 
less than 0.2 % (k= 2). Validation of the T-SIRCUS facility's ability to 
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achieve its estimated uncertainty is established through 
measurements of a primary standard gold-point blackbody by an 
optical pyrometer that had been calibrated in T-SIRCUS with an 
expanded (k = 2) uncertainty less than 0.1 % [14]. 


10. DISCUSSION 


The SNPP VIIRS instrument doesn't have a radiance responsivity 
uncertainty requirement. SNPP VIIRS is a follow-on instrument to the 
MODIS instruments, and it is assumed that the SNPP VIIRS sensor has 
a similar radiance responsivity uncertainty requirement as MODIS. 
For heritage MODIS sensors, the required laboratory calibration 
uncertainty for radiance responsivity was 5 % or less, k= 1. Both 
laboratory calibration approaches (SpMA/SIS-based and T-SIRCUS- 
based) meet that uncertainty requirement. Looking to the future, 
Decadal Survey missions [15] have increasingly stringent TOA sensor 
responsivity uncertainty requirements. For example, the Ecosystems 
instrument on the Pre-Aerosol-Cloud-Ecosystems (PACE) mission, 
tasked with measuring organic material in near-surface ocean layers, 
has a laboratory calibration requirement of 2 % [16]. The Climate 
Absolute Radiance and Refractivity Observatory (CLARREO) [17] isa 
climate-focused mission calibrated and planned to operate in 
reflectance mode; its foundation is the ability to produce highly 
accurate and trusted climate data records. The reflectance accuracy 
required to fulfill its mission requirements is 0.3 %, k= 2, in the 
visible /near-infrared spectral region. 

The uncertainty for ocean color measurements is 0.5 % and on- 
orbit, vicarious calibration techniques have historically be used to 
achieve the required uncertainties. For example, the standard 
deviation of the sorted Sea-Viewing Wide Field-of-View Sensor 
(SeaWiFS) ocean color vicarious calibration data set (6 years of data), 
using only the subset of "good" data and employing the Marine Optical 
Buoy (MOBY) vicarious calibration data set coupled with sensor 
response trending using the United States Geological Survey (USGS) 
Lunar Model, is approximately 1 %. Extrapolating from the SeaWiFS 
example, it takes between 20 and 40 "good" vicarious calibration data 
points for the mean vicarious gain to converge to within 0.1 % of the 
mission-averaged value [18]. Given monthly lunar views, it takes a 
minimum of 2 years of data for the MOBY-based vicarious calibration 
to converge to within 0.1 % of its mean value. Reprocessing of the data 
set is required to apply the vicarious calibration-derived gain 
coefficients, unless the first 2 years of the mission data set is discarded. 
Vicarious calibration does not fit within the current VIIRS operational 
paradigm, which does not allow for reprocessing of data. 
Consequently, an alternate on-orbit approach, without the multi-year 
vicarious calibration latency and potential loss of data, is desirable. 

Because the best estimate for the radiance uncertainty in lamp- 
illuminated integrating spheres is 3 %, k= 1, in the visible and near- 
infrared regions [13], there is no clear path forward using 
conventional lamp-based integrating sphere sources with current 
calibration and validation approaches to achieve the laboratory 
calibration uncertainty required by the PACE, CLARREO, and other 
future Decadal Survey missions. Potential reductions in the 
uncertainty in the SIS radiance could be achieved using transfer 
standard spectrographs [19] moving from a source-based radiance 
scale to a detector-based radiance scale. The expected reduction in 
the uncertainty in the source radiance is an order of magnitude, from 

= 3%to ~ 0.3 % in the Vis/NIR spectral range. 

Using T-SIRCUS to calibrate instruments with uncertainties that 
meet future laboratory requirements is a second option. The T- 
SIRCUS-based laboratory calibration presented here has an expanded 
uncertainty of 0.56 % (k= 2) for SNPP VIIRS with an estimated 
achievable uncertainty of 0.2 % (k= 2). The numbers are notable 
because they represent the first laboratory calibration of an 


instrument that can meet the 0.5 % (k=1) on-orbit ocean color 
instrument uncertainty requirements. 

Whether it is worth the effort to achieve those uncertainties in the 
laboratory calibration of satellite sensors depends on the uncertainty 
in the transfer of the radiometric scale to orbit, the ability to trend the 
sensor response on-orbit, and the ability to account and correct for 
OOB biases arising from the spectral distributions in Earth scenes. 
The following 3 sections deal with the transfer-to-orbit uncertainty, 
effects of source spectral shape on the OOB contribution to a 
measurement, and an approach to correct the OOB signal using a 
variable user-defined virtual calibration source with the T-SIRCUS- 
derived band ASRs. This conceptual approach to the vicarious 
calibration of ocean color sensors using a low uncertainty laboratory 
calibration combined with a low uncertainty on-orbit lunar calibration 
is explored as a complementary or validation approach to the 
vicarious calibration of ocean color sensors using MOBY. 


10.1. Transfer-to-Orbit and Sensor Temporal Changes in 
Responsivity 

An initial attempt to determine the transfer of a satellite’s 
laboratory calibration to orbit was performed by the SeaWiFS project 
[20]. The effect of the transfer-to-orbit was estimated by measuring 
the instrument outputs while viewing a solar diffuser illuminated by 
the Sun at the instrument manufacturer's facility and again at the start 
of on-orbit operations. The ground-based measurements were used 
to predict the instrument's outputs on-orbit with an estimate 
combined standard uncertainty of 3 %. For the 8 bands, the results, 
extrapolated to Day 1, showed the measured outputs averaged 0.8 % 
greater than the predicted signals, with a standard deviation of 0.9 % 
(k=1). The maximum difference between the measured and predicted 
on-orbit signals was 2.1 % for SeaWiFS band 3; the minimum 
difference was -0.9 % for band 1. The uncertainty in the SeaWiFS 
transfer-to-orbit experiment is low enough to support current satellite 
sensors, including MODIS and VIIRS, but is much greater than 
required for future requirements. Clearly, an alternate strategy needs 
to be developed for uncertainties in the transfer-to-orbit 
commensurate with future on-orbit uncertainty requirements to 
justify a requirement for the lowest possible T-SIRCUS-based 
laboratory calibration uncertainty. 

The USGS Lunar Model is based on 6 years of ground-based 
measurements of the Moon from Flagstaff, AZ. Combining monthly 
measurements of the Moon with the USGS Lunar Model, the SeaWiFS 
Project demonstrated the ability to measure sensor band temporal 
changes at the 0.1 % level [18, 21]. While the USGS Lunar Model 
works extremely well as a relative source for trending, estimated 
uncertainties in the absolute lunar irradiance are between 5 % and 
10 % in the Vis/NIR spectral region, too large to support a low 
uncertainty laboratory calibration requirement if lunar 
measurements were to be used in a low uncertainty transfer-to-orbit 
approach. 

Recent hyperspectral lunar irradiance measurements from Mount 
Hopkins, AZ have an estimated uncertainty in the absolute lunar 
irradiance at the top of the atmosphere of 0.5 % over most of the 
instrument's spectral range, from 500 nm to 1000 nm [22]. An 
examination of the uncertainty budget reveals that the dominant 
uncertainty component is the telescope/spectrograph calibration 
uncertainty. Recent measurements imply that moving to a detector- 
based spectrograph calibration could reduce the spectrograph 
calibration uncertainty component from approximately 0.4 % to 
approximately 0.1 % (k= 1) [19]. 

Changing the lunar measurement site location from Flagstaff, AZ, to 
either a high-altitude platform such as an aircraft or a mountain site 
with atmospheric characterization capabilities like NOAA’s Mauna Loa 
Observatory, HI, could potentially reduce the uncertainty in ozone and 
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aerosol concentrations. Ignoring biases due to sensor channel out-of- 
band response (see Section 10.2), allowing for a (k= 1) 0.1% 
uncertainty in band responsivity, a 0.1 % uncertainty in the updated 
USGS Lunar Model, an additional 0.2 % uncertainty in the transfer-to- 
orbit and a 0.1 % uncertainty in the temporal trending of the 
instrument response, the combined standard uncertainty in an on- 
orbit system response is approximately 0.25 %. If achieved, the low 
uncertainty in the absolute lunar irradiance fosters a potential new 
paradigm for establishing the uncertainty in the transfer-to-orbit, 
namely measuring the absolute lunar irradiance on-orbit. 


10.2. Effects of Source Spectral Shape 


For the RSR of Band M11 (see Fig. 14a), there is a strong IB response 
peak around 410 nm and an additional weak OOB response between 
700 nm and 1000 nm. The OOB contribution to a measured 
broadband signal reduces the spectral purity of the band and is 
treated as a measurement error. In practice, there is both a calibration 
spectral distribution and a measurement spectral distribution in the 
measurement equation — whether the instrument is used as a 
radiometer or a reflectometer. The OOB measurement error is 
proportional to the spectral difference between the calibration source 
and the measurement source weighted by the OOB instrument 
response. Ifthe calibration source and the measured source have the 
same spectral distributions, the OOB measurement error cancels, 
independent of the magnitude of the sensor band’s OOB response. 
Consequently, to minimize on-orbit measurement errors due to OOB 
response, it is desirable to have the calibration source spectral 
distribution approximate the measured source spectral distribution, 
to the extent the source measurement distribution is known. 

For the VIIRS bands in the visible and near infrared, the 
fundamental on-orbit measurement is reflectance, which is the ratio of 
the up-welling Earth radiance to the radiance from an on-board 
diffuser irradiated by the Sun. For measurements of radiance, the 
fundamental measurement is the ratio of the up-welling Earth 
radiance to the radiance of the calibration source. Instruments are 
calibrated in the laboratory for spectral radiance responsivity pre- 
launch against lamp-illuminated integrating sphere sources that have 
spectral distributions approximating a blackbody distribution within a 
range of temperatures from approximately 2700 K to 3100 K. The 
SIS100 radiance spectrum used in the historical radiance responsivity 
calibration of VIIRS is shown in Fig. 21a. Figure 21b shows the solar 
irradiance spectrum from the measurements of Thuillier et al. [23] 
used in VIIRS reflectance measurements. Figure 22 shows the IB to 
TB ratios calculated for VIIRS bands M1 through M7 for the two 
spectra shown in Fig. 21. The SoMA-based RSR was used for the IB 
and TB calculations. The effect of the red leak in the band M1 OOB 
combined with the red-rich SIS spectrum create a ratio of 0.76 for the 
410 nm band in (Fig, 22). 
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Figure 21. Comparison ofa laboratory and an on-orbit source spectrum. 
a) SIS100 radiance spectrum used in the laboratory to calibrate SNPP 
VIIRS band M1. b) Solar irradiance spectrum from the Thuillier et al. 
[19] solar model. Note the marked difference in their respective 
spectral distributions. 


Depending on the approach, either the solar spectral irradiance or 
the lamp-illuminated SIS radiance may be the reference source 
spectral distribution. Neither of these spectral distributions 
approximates the spectral distributions of Earth scenes viewed on- 
orbit. Extending the discussion to consideration of Earth scenes, 
Figure 23 shows nominal on-orbit at-sensor radiance distributions for 
a blue ocean scene, a grassland scene, and a desert scene with and 
without trace gas absorption [24]. The symbols give the radiances 
without effects of atmospheric trace gas absorption; the dashed lines 
include atmospheric absorption based on the nominal atmospheric 
absorption described in the reference. They provide reasonable first- 
order estimates for the range of at-sensor spectral radiance 
distributions to be measured on-orbit by SNPP VIIRS. 

Figure 24 shows the ratios of the IB to TB responses for the three 
Earth scenes in Fig. 21 relative to those for a solar calibration 
spectrum and a 2856 K blackbody spectrum. SpMA-based RSR’s were 
again used for the calculations. Results in Fig. 22 are presented as the 
ratio of the IB to TB measurement divided by the ratio of the IB to TB 
calibration: 


| exe | /| DN(IB) | 
Ordinate = | ————~ meas ————~ cal 
DN(TB) DN(TB) . (14) 
where meas is the on-orbit measured spectrum (from Fig. 23) and cal 
is the calibration spectrum (from Fig. 21). Here the measured spectra 
from Fig. 23 are without trace gas absorption. Figure 24 shows the 
errors due to OOB response from using calibration sources with 
spectral distributions that are different from the at-sensor radiance 
spectra on-orbit. In general, using a solar calibration source 
distribution creates residual out-of-band errors equal to or smaller 
than a blackbody calibration distribution. The error from the red leak 
in the band M1 OOB combined with the red-rich blackbody source 
spectrum is particularly prominent in Fig. 24a. However, using the 
solar diffuser as the on-orbit calibration source, OOB errors for the 3 
scenes examined are less than the reflectance calibration requirement 
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Figure 22. The In-Band to Total-Band ratios (IB/TB) for bands M1 
through M7 for a source with a solar distribution and a source witha 
2856 K distribution. These are the errors arising from the OOB 
responses in the bands. The source spectra are shown in Fig. 1. 


The VIIRS sensor bands have been chosen to view Earth-leaving 
radiances in atmospheric absorption windows, that is, in spectral 
regions where there is little absorption from trace gases. As a result, 
trace gas absorption mainly impacts the OOB measurement and 
uncertainties in the at-sensor spectral distribution from these 
absorptions are relatively small. Figure 25 shows the ratio in the 
calculated OOB error, both with and without absorbing gases in the 
atmospheric path. The additional error from trace gas absorption is as 
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Figure 23. Three nominal Top-of-the-Atmosphere (TOA) radiance 
spectra [24]. The symbols give the model TOA radiances without trace 
gas absorption. The dashed lines give the radiances with trace gas 
absorption included. See the reference for details on the derivation of the 
a) spectra.TOA radiances over a blue ocean; b) TOA radiances over a 
grassland; and c) TOA radiances over a desert. 
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Figure 24. IB/TB ratios for three Earth scenes relative to a solar 
calibration spectrum (open symbols) and a 2856 K blackbody spectrum 
(closed symbols). The ratios show the errors arising from the use of 
calibration spectra that are different from the spectra measured on orbit. 
a) Full Scale.; b) Expanded scale. 


large as 0.7 % for measurements of grasslands at 410 nm. The 
maximum error for a desert scene is 0.4 %, also for band M1. The blue 
ocean scene is not sensitive to the atmospheric absorption in this 
example, with additional errors less than 0.1 % for all 7 bands. 
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Figure 25. Effect of trace gas absorption on the measurements (ratio of 


with absorption to without). These are the errors due to trace gas 
absorption. 


Because of the large biases that may be introduced ina 
measurement due to different calibration and measurement source 
spectral distributions, most Earth remote sensing multi-band filter 
radiometers, e.g. MODIS or VIIRS, have OOB corrections algorithms. 
Reference [26] describes one approach applied to SNPP VIIRS. In 
Section 10.3 we presenta slightly different approached based on the 
absolute calibration of the VIIRS filter bands, using the conventional 
measurement approach dividing a measurement spectrum by a 
calibration spectrum. 


10.3 Source Spectral Shape Ocean Color Example 


Science requirements are strictest for measurements of ocean 
color, requiring an on-orbit radiance responsivity uncertainty of 0.5 %. 
One of the advantages of the T-SIRCUS ASR calibration is that any 
source distribution can be used as a proxy calibration source in Eq. 12 


to develop a band-averaged radiance L;" (SiRCUS ) and band-averaged 
calibration coefficients R,, (SIRCUS ) from Eq. 7. This is a user-defined 


virtual calibration source spectral distribution 7°“ (7) ; there is 


‘user —defined 
no uncertainty in the distribution. DNs are calculated from the 
summed product of the ASR and the user-defined calibration source 
distribution. 

To evaluate the impact of finite OOB response from VIIRS bands, a 
matrix of at-sensor radiances (nadir view) was developed using 
MODTRAN [25] for a range of column water vapor concentrations 
(0 cm to 5 cm) and near-surface phytoplankton chlorophyll-a 
concentrations (0.05 mg/L to 3 mg/L) covering the range of values 
expected on-orbit. A standard marine aerosol, ozone at 300 Dobson 
Unit (DU), and a standard pressure were selected and not varied 
during the simulations. A set of 275 at-sensor radiances with differing 
atmospheric and oceanic radiometric properties was calculated for 
SNPP VIIRS focal plane detector 8 and bands M1 through M7. 

Results for 3 different calibration source spectral distributions were 
compared with ‘true’ radiances and histograms of out-of-band biases 
were generated. In addition to the typical calibration sources, an 
incandescent lamp operating at 2856 K and a solar irradiance 
spectrum [19], a proxy calibration source spectrum consisting of the 
modeled at-sensor radiance for a water-leaving radiance spectrum 
with a near-surface phytoplankton chlorophyll-a concentration of 0.5 
mg/Land a column water vapor concentration of 1 cm was used in 
the simulations. This is the virtual calibration source distribution to be 
used with the T-SIRCUS band ASRs in deriving T-SIRCUS-based 
calibration coefficients. 
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Results of the simulations are shown in Fig. 26; the biases arising 
from OOB response are clearly shown for each calibration source. 
Each band shows a static bias and a width resulting from the different 
modeled at-sensor radiance distributions. It is noteworthy that using 
a nominal ocean/atmosphere calibration source spectral distribution, 
biases in VIIRS bands are essentially reduced to 0, with the 
distribution widths of all channels negligible with the exception of 
Band M4, which has a width < + 0.05 %. 

As demonstrated by Gao and Chen [26], considering the suite of 
bands and including any information about atmospheric absorption 
will reduce these residual errors further. Counting the medium 
resolution and imaging bands, there are a total of 10 visible and near- 
infrared bands on VIIRS and some a priori knowledge of the spectral 
distributions of Earth- reflected radiation. This information can be 
used to develop a nominal Earth at-sensor spectral distribution for 
each pixel and from there calculate band responsivities for that 
distribution. With even a simplistic algorithm, allowing flexibility in 
the distribution ofa calibration source, the residual OOB errors may 
be reduced to a negligible level where OOB biases would not impact 
the T-SIRCUS-based radiometric uncertainties propagated to orbit. 


11. CONCLUDING REMARKS 


The traditional broadband source-based calibration of SNPP VIIRS 
is compared with the laser-based T-SIRCUS calibration of the sensor. 
For the radiance calibration coefficients (responsivities), the 
traditional coefficients average 2.2 % higher than T-SIRCUS (+ 1.7 %, 
k=1). The estimated uncertainty for the traditional method is 3 % 
while the uncertainty budget for T-SIRCUS is 0.56 % (k= 2); and the 
calibration coefficients agree within the combined uncertainties for 
the two methods. For the center wavelengths of the central peaks of 
the bands, the traditional values average 0.06 nm lower than T- 
SIRCUS (+ 0.14 nm, k = 1), again within the combined uncertainties for 
the two methods. 

The radiance uncertainty for the T-SIRCUS measurements of SNPP 
VIIRS is more than a factor 10 improved over the traditional method. 
In addition, the 3 % uncertainty for the traditionally measured 
responsivities represents the state-of-the-art for that method. An 
uncertainty analysis has shown the potential to reduce the T-SIRCUS- 
based expanded (k = 2) calibration uncertainty from 0.56 % to less 
than 0.2 %. The T-SIRCUS wavelength uncertainty is more than a 
factor of 5 improved over the traditional method. 

The VIIRS instrument experienced crosstalk between the bands. 
Since the spectral response measurements for the traditional method 
illuminate a single band at a time, crosstalk in the derived RSRs is 
minimal and an extensive set of additional measurements were 
required to estimate the magnitude of this crosstalk. The T-SIRCUS 
technique floods the entire VIIRS focal plane for each laser wavelength 
as the Earth radiance scenes do on-orbit and crosstalk is evident in the 
derived RSRs. The T-SIRCUS measurements of crosstalk confirmed 
that the magnitude of the artifact on an on-orbit measurement is 
small; the results have relieved concerns of the ocean color 
community over the potential impact of large optical cross-talk on 
SNPP VIIRS data products 

With respect to determining detector-to-detector responsivity 
differences within individual bands, the measurement of all detectors 
simultaneously by T-SIRCUS provides an improved characterization 
over the traditional method. From a preliminary examination, the T- 
SIRCUS-determined detector-to-detector differences agree well with 
the effects seen on orbit. 

The T-SIRCUS calibration suffered from lack of flux in several 
spectral windows, limiting the dynamic range of OOB measurements. 
So, while the VIIRS responsivity is measured with low uncertainty in 
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Figure 26. Histograms of out-of-band measurement biases arising from 
different VIIRS calibration sources, shown in bold type on the top left 
hand side of each graph. A total of 275 spectra with different 
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atmospheric and oceanic properties were used to develop the histograms. 


spectral windows where the T-SIRCUS lasers output sufficient power 
to give the required OOB dynamic range, the T-SIRCUS measurements 
could not fully replace the OOB SpMA measurements. In the interim 
since the SNPP VIIRS calibration, a new laser system has been 
developed that rectifies issues with the T-SIRCUS system, including 
the low radiant flux in the 500 nm to 565 nm spectral region. The new 
system was used in the radiometric characterization and calibration of 
the JPSS J1 VIIRS sensor in December 2014-January 2015 [27]. 

The reduced uncertainties in the laboratory calibration of Earth- 
observing satellite instruments using T-SIRCUS will require additional 
improvements to on-orbit uncertainties to be effective. These include 
a reduction in the uncertainty of the transfer of the laboratory 
calibration to orbit as well as a minimization of errors due to OOB 
response for varying Earth scenes on orbit. Suggestions for both of 
these effects were discussed. The characterization of long-term 
instrument responsivity changes can be addressed by monthly 
measurements of the Moon [21] while the transfer-to-orbit 
uncertainty component may be based on measurements of the 
absolute lunar irradiance [22]. With the resolution of these effects it 
may be possible to meet the calibration requirements for many future 
instruments listed in the Decadal Survey [15], though meeting the 
CLARREO mission reflected solar uncertainty requirements will 
remain challenging. 
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